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CHAPTER I 
INTRODUCTION 

1.1 Motivation an d Relevance of Report. In April of 1974 
the National Aeronautics and Space Administration conducted 
a high altitude balloon experiment called LACATE (Lower At- 
mosphere Composition and Temperature Experiment) which em- 
ployed an infrared radiometer to sense remotely vertical 
profiles of the concentrations of selected atmospheric trace 
constituents and temperature. The constituents were measured 
by inverting infrared radiance profiles of the earth's hori- 
zon. The radiometer line of sight was scanned vertically 
across the horizon at approximately 0.25° per second, requir- 
ing 30 seconds to acquire a complete radiance profile. The 
specifications require that the relative vertical position of 
the data points making up a profile be known to approximately 
30 arc seconds. 

The balloon system for accomplishing the mission (see 
figure 1) consisted of: (1) a 45 million cubic feet balloon, 

(2) a load-bar containing the balloon control equipment, (3) 
a package containing additional balloon control electronics 
and a gondola recovery parachute, and (4) a gondola containing 
the research payload. Instrumentation to determine altitude 
consisted of a magnetometer and three orthogonally oriented 
precision rate gyros, the latter were used to obtain accurate 
time histories of the roll, pitch and yaw motions of the gon- 
dola. The magnetometer and gyros were flown with the research 
payload and their outputs were telemetered to ground operations 
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for recording and real time data reduction and display. 

Figure 2 gives the dimensions and weight of the various 
subsystems comprising the balloon system. Two kinds of conn- 
ectors between subsystems were employed: a clevis and ring 

shown in figure 3a and a triangular plate connection shown in 
figure 3b. 

In order to fix the orientation of the line of sight of 
the radiometer, it is necessary to be able to determine the 
configuration of the platform in space, i.e. the attitude of 
the system. This can be accomplished by simulating the bal- 
loon system and using the gyro output in conjunction with a 
parameter estimation process. This process is described in 
greater detail in reference l . The critical problem arising 
in the simulation of the high altitude balloon system is the 
development of the mathematical model. This model must enable 
one to predict the orientation of the platform with sufficient 
precision such that the position of the data points can be 
determined within the required values. At the same time the 
model must be as simple as possible in order to be amenable 
to simulation. 

The actual motion of the balloon system (once it reaches 
float altitude) is extremely complex and involves various 
types of oscillations including bounce (vertical oscillation) , 
pendulations (in plane motion) and spin (rotation) . In the 
general case, where the system is subjected to arbitrary in- 
put (initial conditions) these oscillations will be coupled. 
The modeling problem is further complicated by other factors; 
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e.g. the presence of damping forces, the significant mass 
(relative to the payload) and elasticity of the cable, the 
tendency (due to support conditions) of some of the gondolas 
to oscillate independently of the main system, and the fact 
that the balloon itself is actually a distributed parameter 
system and hence its properties (e.g. moment of inertia) 
must first be approximated in order to develop a lumped para- 
meter model for the entire system. 

1.2 Objectives of Report. The main objective of this report 
will be to develop a mathematical model for predicting the 
three dimensional motion of the balloon system. The model 
will incorporate the various factors discussed previously and 
includes the effects of bounce, pendulation and spin of each 
subsystem. Also, some work will be done in analyzing the bal- 
loon itself, i.e. determine the best way for treating it as 

an equivalent lumped parameter system. This will require some 
investigation of boundary layer effects and the aerodynamic 
forces acting on the balloon. Finally, various simplified 
forms of the system mathematical model will be developed based 
on an "order-of-magnitude" analysis. 

1.3 Idealization of System . In general, it is necessary to 
idealize the physical system before one obtains the mathemat- 
ical model. For purposes of this study, the following assump- 
tions will be made: 

(1) The masses of balloon, subsystems and interconnecting 
subsystems will be "lumped" at the locations shown in figure 4. 
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(2) The balloon will be treated as an "equivalent" 
rigid body. 

(3) The generalized drag reactions will be computed by 
assuming that balloon and masses are spherical in shape. 

(4) The altitude of the balloon static equilibrium pos- 
ition (float altitude) will be assumed to be a constant during 
the entire period of observation, i.e. changes in this alti- 
tude due to the losses of or changes in the properties of he- 
lium will be neglected. 

(5) The cables will be treated as though they were in- 
flexible . 

Further assumptions are made in this study but these will 
be discussed in the main body of the thesis. 
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CHAPTER II 
BASIC CONCEPTS 

2.1 Equilibrium Principles of ffechanics . The vector (Newton's law) 
and the variational (Hamilton's principle and Lagrange's Equation) 
principles of mechanics provide two methods for obtaining mathematical 
descriptions of the same realm of natural phenomena. The vector theory 
bases everything on two fundamental vector quantities (momentum and 
force) , while tire variational theory bases everything on two scalar 
quantities (kinetic energy and generalized work function) . In the 
case of free particles, i.e. particles whose motion is not constrained, 
the two forms of description lead to identical equations. However, 
for systems with oons train ts the analytical treatment is simpler and 
more efficient. With the variational method, the given constraints are 
considered in a natural way by letting the system move along all tire 
possible paths (in configuration space) in harmony with them. With 
the vector method, the forces which maintain the constraints must be 
considered. The vector method does not restrict tire nature of a force, 
while the variational method requires that the acting forces can be 
derived from a scalar quantity, termed tire work function. Forces, 
which cannot be derived from a generalized work function (e.g. Coulomb 
friction) must be included via the generalized force concept. 

There are two advantages in modeling tire balloon system under 
study by the analytical method. 

(a) Tire energy terms and generalized forces can be computed in 
a straightforward manner without great difficulty. 
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(b) Hie internal reaction forces which exist at the interface 
of the subsystems do no work during the motion, and hence 
they can be ignored. This results in a simplification of 
the form of the model. 

2.2 Hamilton's Principle and Lagrange's Equation . Hamilton's 
principle states that the motion of any mechanical system in configura- 
tion space is such that 



(T + W*)dt = 0 , 


( 2 - 1 ) 


where: T = kinetic energy, 

n 

W* = E Q.q. = generalized work function, 
i=l 33 

n = number of degrees of freedom, 
n _ 3r . 

. Q. = E F.*— — = generalized force, 

3 i=l 1 dq i 

i 

F. = external forces acting on system, 

r. = position vector to point where the external force acts, 
q. = generalized coordinates, and 
N = number of external forces. 

From the calculus of variation, the necessary and sufficient 
conditions (Euler-Lagrange Equations) for equation (2-1) . to be 
satisfied are 


_d_ 3T _ 9T_ 
dt 9^7 9q ± 


(i 1. . .n) , 


( 2 - 2 ) 


and Sq^tf) = 6q^(t 2 ) =0. 

Moreover, if (a) the generalized forces can be divided into conserva- 
tive and non-conservative components, i.e. 
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0- = Q. C + p. N , where 

l l l ' 

Qj~ = 9V/9q i , and if 

(b) the potential energy (V) is not a function of the generalized 
velocities (q^) , then equation (2-2) can be written as 


d 9L 9L _ N 

dt 9?^ 9q^ ^i ' 


(2-3) 


where the Lagrangian L = T - V. 

Equation (2-3) can be written in state- variable form by first 
defining a generalized momentum (p^) , 



Substituting (2-4) into (2-3) yields 


9L , _ N 

p. = - — + 0. 

1 1 3q t 1 


(2-5) 


Equations (2-4) and (2-5) can be employed to obtain tlie 2n state 
variable equations (in terms of generalized coordinates and generalized 
momentum) by inverting equation (2-4) to get tlie generalized velocities 
in terms of p^ and and then substituting these expressions into 


N 


equation (2-5) . 

If (a) tlie system is conservative (Q.” = 0) , 

(b) tlie generalized forces can lie divided into applied 

- 

forces (Q. ) and reaction forces (Q. ) such that 

i l 

(c) tlie work done by tlie reaction forces during a virtual 
displacement is equal to zero, and 

(d) tlie applied forces can be written in tlie form 
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0 a = .^L 

Q i 3q i ' 


then, equations (2-4) and (2-5) yield the Canonical equations of 
Hamilton; i.e. 


* 




( 2 - 6 ) 




9H 


9 p. 


where H, the Hamiltonian, is defined by 

II(p,q,t) = Eq i P i - L(q i ,q i ,t) . 
i 


(2-7) 


2.3 Generalized Coordinates . The generalized coordinates for a given 
system are those coordinates which are employed to specify the configura- 
tion of the system at any instant of time. In any mechanical system 
there will be as many generalized coordinates as there are degrees of 
freedom. In tlie case of the idealized lunped parameter system shown 
in figure 4, eleven generalized coordinates are required to specify 
tlie configuration. These consist of three translational coordinates 
to locate tlie mass center of the balloon (relative to a set of axes 
fixed in the space) and eight Euler angles to specify tlie orientation 
of tlie three pendulum subsystems. 

In general, the Euler angles give tlie orientation of tlie body 

coordinate system (x ) relative to a fixed system (X.). If tlie two 

i 

systems are initially coincident, a series of tliree rotations about 
the body axis is sufficient to allow tlie tody axes to attain any 


orientation. 
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Several sets of Euler angles are possible for fixing the orientation 
of the balloon and subsystems. TI iree sets were investigated in this 
work and they are shown in figure 6. ihe sequence of the three 
rotations which define the set which was finally employed in this work 
is described be lav. 

(a) A positive rotation 0 about the axis resulting in the 
x| body system, 

(b) a positive rotation ip about the x^ axis resulting in the 

a- 

xh body system, and 

(c) a positive rotation <J> of the tody about the x^ axis relative 
to the xV system. Since the final rotation is relative to 
the xV system, the transformation equation is given as 

x" = A X , (2-8) 

where 

c(<jO s (ip) s ( 0 ) -s(ip)c(e)' 

0 c ( 0 ) s (0) (2-9) 

s(ip) -s(e)c(ifO c(4»)c(e) 

This set is convenient whenever a body has at least one axis of 
symmetry. 

2.4 Aerodynamic Reactions . Hie exact mathematical model for the 
balloon system consists of the equations of motion of the solid 
(balloon fabric and payload) and the fluid-dynamic equations. These ■ 
equations are coupled through the boundary conditions which must be 
satisfied at the interface of the solid and fluid media. The resulting 
mathematical model is extremely corrplex and consists of a system of 
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coupled non-linear differential equations which must be solved 
simultaneously in order to predict the behavior of the system. 

For purposes of this work the mathematical model will be simplified 
by treating the system as a lumped parameter system. In order to 
accomplish this it is first necessary to develop approximate expressions 
for the forces and torques which result due to the interactions between 
the solid system and tie fluid media (air-helium) . In tie lumped 
parameter model, tiese forces and torques are tien treated as external 
reactions on the solid system. 

The reactions which tie fluid media exert on tie solid system 
consist of (a) a "static" lift force due to buoyancy, (b) "dynamic" 
buoyancy forces due to the acceleration of tie fluid approaching tie 
system, (c) viscous drag forces due to tie translational motion of a 
solid body in a fluid medium, (d) viscous 'drag torques due to tie 
rotation of tie balloon in tie fluid media, and (e) inertia drag forces 
and torques due to tie unsteady motion of tie fluid media relative to 
the balloon system. Hie expressions for these reaction forces and 
torques are presented below. 

2.5 Static Buoyancy Force . Hie net lift force due to tie aerostatic 
pressures is given as 

P L “ V lPK * ll’ ' <2 - 10) 

where F = static lift force in tie z direction, 

Jh 

V„ = volume of helium, 
p = density of air, and 
P n = density of helium . 

If we assume (a) that the mission takes place in tie isothermal 
region of the atnpsphere and (b) tiat helium undergoes an adiabatic 
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process, then the expressions for and are given as 


and 



II 



_1_ 

"il 


r 


( 2 - 11 ) 


( 2 - 12 ) 


where p Qa = density of air at the static equilibrium position 
(i.e. the float altitude) , 

p q jj = density of helium at the static equilibrium position, 

P o = atmospheric pressure at static equilibrium position, 
z = elevation measured from the static equilibrium position, 
and 

n^ = ratio of specific heats for helium, (c ^/c^) . 

Substitution of (2-11) and (2-12) into equation (2-10) yields the 
expression for the static lift force (F T ) as 

l_i 


P H V H 9 


oa 


p oa^ z 

Po 


oil 



(2-13) 


The presence of a static-buoyancy force results in an external 
torque (see figure 7) whenever the balloon shape is such that the 
center of gravity of the shell does not coincide with tire center of 
buoyancy. The magnitude of this torque is given as 

N l = d x F l , (2-14) 

where = external torque, and 

d - vector between the center of gravity and center of 
buoyancy. 
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Equation (2-16) indicates that the static lift force is a non- 
linear function of z. A linear expression for F , valid for small z, 

ij 

can be obtained from a truncated Taylor-series expansion: 


F = F + 
L o 


fdF L l 


dz 


z=0 


where F q = lift force at the static equilibrium position. 

<F 0 = V) , 

w = weight of the solid balloon system. 

dF L 

The equivalent spring constant I< = -= — __ is obtained by 

o az z — vj 

differentiating equation (2-13) and considering p^V^g to be a constant. 
This development is presented in Appendix A • 


2.6 Dynamic Buoyancy Force . The absolute acceleration of the fluid 
approaching tire balloon results in a uniform pressure drop across the 
balloon in the direction of the acceleration of tire fluid. This 
pressure drop, when integrated over the balloon surface, results in a 
force which acts on the balloon in the direction of tire absolute 
acceleration of tire fluid. The expression for this force, termed 
dynamic buoyancy force, is given as 

h - vA (2 - i5) 

where F„ = dynamic buoyancy force, and 

V„ = absolute acceleration of air. 

W 

This force acts at tire center of buoyancy and results in a torque 
whenever the center of buoyancy does not coincide with tire center of 
gravity (see figure 7) . 
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2.7 Inertia Drag Force . Whenever a body accelerates in a fluid 
nedium, it is subjected to a resisting force even if the fluid is 
assumed to be frictionless. This resisting force, termed the inertia 
drag force, has the same effect as that obtained by increasing the 
body's mass by an amount equal to that of the fluid carried along with 
it. The expression for the inertia drag force on a spherical body 
due to the external flew of air ( 2 ) is given as 


F 


I 


- - p V„ 

2 a H 



(2-16) 


v/hile, for the internal flow of helium, the expression is given as 


F I = - 


p H V H V b 


¥ 


where = inertia drag force, and 

V. = acceleration of the balloon. 
, b 


(2-17) 


2.8 Viscous Drag Force . The egression for the resisting viscous 
drag force ( 3 ) wloich acts on a spherical solid due to its translation 
relative to air is given as 



( 2 - 10 ) 


where = viscous drag force, 

C D = drag coefficient, 

A = projected area of the body on a plane normal to the 
relative velocity, 

V_ = balloon velocity, and 
V w = air velocity. 

The magnitude of depends on the Reynolds nunfoer. 



21 


2 . 9 Viscous Drag Torque due to the Rotation of a Spherical B al loon 

in Viscous Fluid Media . The expression for die drag torque on 

a balloon rotating about a single axis depends on the range of the Reynolds 

wr“ p 

numbers under consideration. For high Reynolds numbers, Re = — - — — , 
and | |tu| 

where 


The above expression is valid for both oscillatory and non-oscillatory 
motion. 

Equation (2-19) can be used to obtain the drag torque due to both 
internal and external fluids and was developed under the following 
assumptions : 

(1) The boundary layer thickness of the fluid is a constant and 
is taken to be equal to the boundary layer tliickness on an 
infinite rotating disk. This thickness (5 ) is given as 


>> | | io | | , the expression ( 4 ) is given as 


T , , . 1/2 1/2 4 *, 

N v = -I<^ (py) a) r o a) 


(2-19) 


Ny = viscous drag torque, 
y = viscosity of the fluid, 

a) = angular velocity of balloon along the rotation axis, 

r = mean radius of the balloon, 
o 

w = angular acceleration of the balloon, 

| co [ | = amplitude of w, and 
Jwj |= amplitude of w. 


5 


4.5 


y 

pe 


r 


*Refer to .Appendix B for development of this expression. 


( 2 - 20 ) 
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where 

6 = boundary layer thickness, and 

co = nagnitude of angular velocity of the balloon. 

(2) The boun dar y layer thickness is assumed to be small compared 
to the mean radius of the balloon, and 

(3) Hie tangential velocity is assumed to vary parabolically 
across the boundary layer thickness. 

For the case of a sphere undergoing oscillatory motion with 


frequency 0 about a single axis and having low Reynolds nunber, the 
expression ( 6 ) for the viscous torque due to the relative motion of 
the inner fluid is given as 

= ' I 17 y H r o 3 (B K r o ) “ ' • C 

where Ny = internal viscous torque due to helium, 
y = viscosity of helium, 

Q = frequency of rotational oscillation, and 


v IT = kinenatic viscosity of helium. 

H 

For the relative motion of the outer fluid, the express ion ( 6 ) for 
the torque is given as 




Try r 
3 a o 


3 + 63 r + 63_ 2 r 2 + 23 V 3 
a o a o a o 


1 + 23 r + 23 2 r 2 
a o a o 


0) 


( 2 - 22 ) 


where N = external viscous torque due to air, 

= viscosity of air, 

3 = (n/2v ) 1//2 , and 

a & 

v = kinematic viscosity of air. 
a 
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2 . 10 Inertia Drag Torque due to the Unsteady Rotation of a Spherical 
Balloon in Viscous Fluid Media . The expression for the inertia 
drag torque on a spherical balloon rotating about a single axis depends 
on the magnitude of the Reynolds number. For high Reynolds number and 
[ | co | | >> (| u [|, the expression is (4) 

N I = " K 2 *' (2-23) 

where = inertia drag torque. 

Equation (2-23) is valid only for the case of non-oscillatory motion 
but can be used to determine the inertia drag torque due to both inner 
and outer flews. 

For low Reynolds nurrber, tire expression ( 6 ) for the inner flow is 
given as 

N I = - I * P H r o 5 (6 H r o ) ~ 1 W ' (2_24) 

where N_ = internal inertia drag torque. 

For the case of outer flow, the expression for the torque is given as 


.. 8 5 

N T = - up r 
I 3 a o 


1 + 

3 . O 


1 ’+ 28 r + 28 2 r 2 
a o a o 


w , 


where N_ = external inertia drag torque. 


(2-25) 


*See Appendix C for development of this expression. 
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chsptEr iii 

DEVELOPMENT OF STATE VAFJAELE EQUATIONS 
3.1 Kinematics . In order to obtain the kinetic energy of the balloon 
system, it is first necessary to develop the kinematic expressions for 
the velocities (angular and linear) of the balloon and subsystems shewn 
in figure 4 . By enploying the Eulerian Angles which were described in 
Chapter II, the angular velocities for the balloon and subsystems can 
be written as 

= u 1 = 0 1 c(^ 1 ) e* + ip 1 e* + (e-jSf^) + e* , 

oj = © 20 (^ 2 ) + ^2 e 2 + (02 s ^ 2 } + ^1^ e 3 ' (3-1) 

—I ' —3 * —3 ' ’ —3 

to = e^c (^3) + ^3 e 2 + ( 0 3 s (^3^ + ^3) e 3 > 

where uP = angular velocity of the balloon, 

th 

u 1 = angular velocity of the i subsystem (i=l,2,3), 

6 . ,4^ ,<Jk = Euler angles of rotation, and 

e? = unit vectors for the subsystem (see figure 4 ) 

For purpose of this work, the spin of subsystem 2 is assumed to be 
identical to the spin of the balloon, i.e. cf>^ = <J> 2 * This assumption is 
reasonable because of the type of connectors employed (refer figure 3 ) • 
The linear motion of the balloon center is referred to an axis 
which is fixed in space at the balloon static equilibrium position 
(see figure 8 ) . The velocity expression is given as 

v® = x i + y j + z k , (3-2) 

where v^ = velocity of balloon center, and 

* ■ • 

x,y,z = components of the velocity of the balloon center. 
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The velocity 'expressions for subsystems are given as 


— i — i-1 — i — i 

v = v - r^ w x e^ , 


where v 1 = velocity of the mass center of the i^ 1 subsystem, 


(i = 1,2,3) 


(3-3) 


th 


^ = velocity of the ness center of the (i-1) subsystem, 

r . = distance between m. and m. , . 

i l l-l 

The detailed velocity expressions in terms of the Euler angles and 
their rates arc obtained by substituting equation (3-1) and (3-2) into 
(3-3) and these are presented in Appendix D ■ 

3.2 Kinetic Energy . The general kinetic energy expression for the 
balloon shell which is undergoing both rotational and translational 
motion is 


m 

o 


V°.V° 


ol 


( w °) 2 + ( co °) 2 


o3 , o, 2 
+ -j- to 3 > 


(3-4) 


where 


m 


= kinetic energy, 

= mass of balloon shell. 


I ,l „ = moments of inertia of tlie balloon shell relative 
ol o3 

to the e? modified body axis, and 
= conponents of the balloon angular velocity along 


eT modified body axis. 

The detailed expression for the kinetic energy is obtained by sub- 
stituting equations (3-1) and (3-2) into equation (3-4). This 
expression is presented in Appendix H • 

The kinetic energy egressions for the subsystems arc obtained by 
substituting equations (3-1) and (3-3) into the following 


i m i _i _i ^1 , i % 2 I i2 , i. 2 , I i3 , i.2 

T = -j- v *v + — (u^) + ~ 2 ~ + ~T (u> 3 ) ' 

(i=l,2,3) 


(3-5) 
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where = kinetic energy of the i^ 1 subsystem, 

nu = mss of the i^ 1 subsystem, 


I. . = moments of inertia of the i^ 1 subsystem relative to 
il 


et modified body axis, and 


til 

uk = components of angular velocity of the i subsystem 
along e^ direction. 

For purposes of this work, the moments of inertia I^and I^were 
neglected in which case equation (3-5) becomes 


m. 


T 1 = ^ V 1 -^ + ( u i) 2 , (i=l,2,3) 


(3-6) 


The detailed egressions for T 1 are given in Appendix E . Tire total 

kinetic energy (T) of the balloon system is obtained by sunming up the 

expressions in equations (3-4) and (3-6), i.e. 

3 
Z 

i=l 


T = T° + Z T 1 


(3-7) 


3.3 Potential Energy . The expressions for potential energy due to the 
gravitational forces acting on the subsystems are given as 


V 1 = - m. g I r .c(i|».)c(0 .) , i=l,2,3 


1 ^ : 


(3-8) 


where V 1 = potential energy of the i Ul subsystem with datum 
through the balloon mss center. 

Tlie expression for tlie potential energy due to the static buoyancy 
force is given as 


' K 0 

V = -S- z 2 
s 2 


(3-9) 
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where V s = potential energy due to static buoyancy force, 

K q = equivalent spring constant, and 

z = elevation measured from balloon static equilibrium 
position. 

The total potential energy (V) for the balloon system is obtained by 
sumtdng up equations (3-8) and (3-9), i.e. 

3 i 

V = X v 1 + V (3-10) 

i-1 

The Lagrangian (L) , which is defined in equation (2-3 ) , is 
obtained by subtracting the total potential energy in equation (3-10) 
from the total kinetic energy in equation (3-7) and is given as 

3 . 3 . 

L = T- V=T°+ XT 1 - I r - V . (3-11) 

i=l i=l 

3.4 Generalized Forces . The generalized force (h, corresponding to 
the generalized coordinate q^, can be determined from equation (2-1 ) 
or, in the case of torques, by determining the virtual work done by the 
torque during a virtual displacement 6q^ (7 ) • In this manner the 
generalized forces corresponding to each of the aerodynamic reactions 
(refer Chapter II) were determined and the resulting non- zero 
expressions are presented below. 

(a) Static Buoyancy Force* 

QqI = w d s(G^) , and (3-12) 

Qjl = w d s(i)> 1 ) . 

*The generalized forces presented here are due (only) to the torque of 
the static buoyancy force, since the force itself is considered via 
the potential energy expression (see equation 3-9 ) . 
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(b) Dynamic Buoyancy Force* 

Q x = vA ' 

<4 = Va’i ’ an3 

Q l m Va’3 ' 


where 


V__ = v.T + v„T + vJc 


(3-13) 


= wind velocity. 

(c) Inertia Drag F orce 

- -0.5p a V H (x - v x ) - c„V H x , 

= -0.5p a V H (y - v 2 ) - p r V h y , and (3-14) 

Q z = -0.5p_V H (z - v 3 ) - p R V H z . 

(d) Viscous Drag Force* * 

Q 4 = -§ CpAp^U-v^ 2 + (y-v 2 ) 2 + (z-v 3 ) 2 ] 1/2 (i-v 1 ) , 

qJ = - i- CpAp^^-v-j^) 2 + (y-v 2 ) 2 + (z-v 3 ) 2 D 1/,2 (y-v 2 ) , & (3-15) 

Q 4 = - | < V^ a [I(x-v 1 ) 2 + (y-v 2 ) 2 + (z-v 3 ) 2 D 1/2 (z-v 3 ) , 

where A = tt r Q 2 . 

(e) Viscous Drag Torque 

(1) High Reynolds Nurrbers 

q 5 q1 = -3.72 r Q 4 C (P a v a ) 1/2 + (p h p h ) 1/2 ^ (®i +| ( , i + ^i) 1/4 (9i + 4>i s (%) ) • 

= -3.72 r o 4 C(P a P a ) 1/2 +(P H %) V2 ](0^^i) 1/4 ^i)^ & (3_16) 

Q®! = -3.72 r G 4 ^ ( P a P a ) 1/2 + ( P H %) 1/2 ^ ( 5 i + ^l + ^ i ) 1/4 < 0 X s ^ i ) + V • 


* The torque due to the dynamic buoyancy force was neglected in this 
study. 

**The translational drag forces acting on the subsystems itu are 
neglected in this study. 
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(2) Low feynolds limbers 


— 

2 2 3 3 


I 

3+68 r +63 r +28 r 

8 3 , a * 

3 "Vo ( 0 hV 

8 3 1 

“ ^ TT]i r ! 

3 a o 

a o a o a o 

1+28 r +28 2 r 2 


^ 3. O 3 O 

— 


( 0 JL + -j^ s ( ) r 



8 3 

- -=■ Tip r 
3 a o 

2 2 3 3] 

3+63 r +63 r +28^r 
a o a o a o 

- § ’Vo 3 Vo 1 



1+28 r +28 2 r 2 

a o a o j 




($>]_) , and 

2 2 „ 3 3l 

C 

Q ^1 = 

8 3 

3+68 r +63 r +23 V 
a 0 a 0 a 0 

8 3 to s 

- 3 "Vo ( Vo } 

“ sr- TTU T 

3 a o 

1+28 r+28 2 r 2 



a 0 a 0 

— 


( 0 ^S (lf^) + <j>^) . 


(f) Inertia Drag Torque 
(1) to >> SI 

Q 0 l = - 6 . 28 C(p a y a )' 1 +(p H v H ) !l ]( 0 1 2 +^ 1 2 +<f> 1 2 ) 



r (6^+^^([i^c(^^)+<()^s (^) ) , 

-6.28[(p a M a ) 2 +(p H u n ) 2 H(0- L 2 +^ 1 2 +<f ) 1 2 ) ' 


r 0 4 ( *l”®l*l c( *l )) ' ^ 

0*1 - - 6 . 28 C(P a li a ) !s +(P Ill . H ) ,S ]( 6 1 2 +i 1 2 4 1 2 )'’ ! 

r o 4 (b l s( h )+ ®l^l c( ^ )+ 4 , l ) " 


i-17) 


1-18) 
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(2) fi >> a) 


<& = 


•r up r 
3 a o 


1+6 r 
a o 


1+26 r +26 2 r 2 
a o a o 


I "Vo 5 'Vo *" 1 


(0^+lf)^4)^C(^)+^S(l|)^) ) r 




8 5 

1+6 r 
a o 

- 5 'Vo 5 ' Vo 1 " 1 

— zr 7TD V 

3 a o 

1+28 r +26 2 r _ 2 | 


& 

o 

0 

— 


) , and 


p *i = 


t- irp r 
3 K a o 


1+6 r 
a o 


1+23 r +26 2 r 2 
a o a o 


3 1Tp Ii r o ($ Il r o ) 


(6 1 s(t/j 1 )+e 1 ij) 1 c(i(; 1 )+<fi 1 ) - 


(3-19) 
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3.5 Generalized Momenta . The generalized momenta corresponding to 
each generalized coordinate is given below 


p =i£ 
x 9x ' 


p = 3T 

y 9y ' 


P = 9T 

Z 9z ' 


P = *£ 

01 


p = *£ 
02 r 


P 

e3 a^ 3 ' 


(3-20) 


p = il 

ij>l 9^, ' 


p = H 

ip2 9 iJj 2 ' 


P.^= 


9T 


i^3 ' 


9T , 

p <j>r 9^' 311(3 


p.,= 


9T 


(f) 3 

The expressions for the generalized nonenta are obtained by substituting 
the expression of kinetic energy in equation (3-7) into equation (3-20) . 
The resulting detailed expressions are presented in Appendix 



33 


3.6 State Variable Equations . The state variable equations are 
obtained by first inverting equation (3-20) and solving for the 
generalized velocity q^ and employing these equations in conjunction 
with the following 


p x - f + JjQx ' 


3L . (j „i 


P y ■ 


K - £ ♦ ■ 


p = 3L + JO 1 

P 0l 39, . , Q 01 ' 

1 1=1 


P 02 = Hr + .=.4 ' 

2 1=1 


p a >L + jo 1 
03 00 3 i=l 03 ' 


• 3L ^ i 

P *1 . - 3^ + 4 Q *1 ' 


P„ = It- + E q), - 

i|)2 9^2 i_i 


p = 1 

>3 9^ 3 + ^*3 ' 


(3-21) 


P., = 


3L + l , and 


*1 " 3+ x 
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CHAPTER IV 

RESULTS AND CONCLUSIONS 

4.1 Linearized Mathematical Model. The general form of the mathemati- 

cal model can be obtained by direct expansion of the state variable 
equations which are presented in Appendix 0 (i.e. equation 3-20 and 

3-21). It is clear that this expansion would be quite tedious and, 
one expects the form of the resulting math model to be extremely 
complex thus requiring a numerical integration method for solution. 

In this work, the form of the mathematical model was simplified by 
neglecting all non-linear terms, i.e. by dropping those terms which 
involve products of the state variables (e.g. q.q. , q- q- , etc.). This 
approximation is justified in view of the fact that actual observations 
indicate that the displacements and velocities of the balloon systems 
are small, i.e. |q.|, (qj « 1. The result of this is a linear model 
which admits a closed form solution thus simplifying the simulation 
process. Two forms of the linearized model were obtained in this study 
and these are discussed below. 

4.2 First Approximation of Mathematical Model . The first approximation 
of the mathematical model was obtained by neglecting: 

(a) viscous drag forces and torques, 

(b) inertia drag torques, 

• • 

(c) the relative spin of the subsystems, i.e. 4> 1 =4 3 2 =< f > 3 , 

(d) the acceleration of tire wind, and 

(e) all nonlinear terms. 

Under these assumptions equation 3-20 and 3-21 can be written 
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P 1 = B 1 q l ' 


P 2 B 2 q 2 ' 


(4-1) 


P 3 B 3 q 3 


P 1 = C i q l + D l q l ' 


P 2 = C 2 q 2 + °2 q 2 r 311(3 


(4-2) 


where 


P 3 C 3 q 3 + D 3 q 3 ' 


P l = 



<P3 ' 



03 ' 


4>1J ' 


q l = 


q 2 = 


q 3 = 


and the elements of the , Ch and matrices (i— 1,2,3) are defined in 
.Appendix H 

The state variable form of the model is obtained by first inverting 
equation (4-1) , i.e. 


q l ~ B 1 P 1 ' 

q 2 = S 2 P 2 ' 
q 3 = e 3 P 3 , 


(4-3) 
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where 


3. = B. 1 
1 x 


(i=l,2 r 3) 


* 


and the B^ ma trices are defined in Appendix I . Differentiation of 
equation (4-3) and substitution into equation (4-2) yields 


i = h % ■ 


P 2 *2 q 2 


P 3 = *3 q 3 ' 


where 


-1, 


A i = (I - Vi> c i 


(i=l,2,3) 


(4-4) 


(4-5) 


The elements of the A^ matrices are presented in Appendix J . Equations 
(4-3) and (4-4) comprise the state variable form of the math model. 

These equations can be written in second order form by differentiating 
equation (4-3) with respect to time and employing equation (4-4) to 
yield 


^1 °1 q l ' 

q 2 — 0-2 ^2 ^ and 


(4-6) 


where 


q 3 a 3 q 3 ' 


a . = 6 .A. 
l li 


(i=l,2,3) 


(4-7) 


The elements of the matrices cu are defined in Appendix K. 

4.3 Second Approximation of Mathematical Model . Hie second approxima- 
tion of the mathematical model was obtained by retaining the dynamic 
buoyancy, viscous drag and inertia drag reactions along with the 
effects of wind velocity and acceleration. However, as in the case of 
the first approximation, all non-linear terms were dropped and the 
relative spin of the subsystems was neglected under these assumptions, 
the state variable equations are given as 
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, and (4-8) 


(4-9) 

as in equation (4-1) , arid 5^ are 
matrices defined in Appendix L and is defined in Appendix M. 

4.4 Eecoirmendation for Future Study . Before enploying the mathematical 
models (refer sections 4.2 and 4.3) in the attitude determination 
process, the models should be evaluated in order to determine tire 
contribution of various terms on the system response and attitude. 

For example, the importance of retaining the inertia drag torque in 
the system model can be determined by: (a) solving tire mathematical 

model given in equations (4-3) and (4-4) , (b) modifying tire form of 
tliis model to include the effect of (see equations (2-23 - 2-25) ) , 
and (c) resolving the modified model. By comparing tire results, e.g. 
natural frequencies, angular velocity, etc., one can determine whether 
this term is important. In tire same manner one can determine whether 
the viscous drag reactions , wind velocity and wind acceleration are 
inportant. Finally tire effect of the non-linear terms can be investi- 
gated by developing a third approximation of tire mathematical model 
which retains first and second order terms and enploying tire same 
process. This work will be conducted in the near future. 


P 2 Y 2 q 2 6+ 2 P 2 + F 2 
P 3 = Y 3 q 3 6+ 3 P 3 + P 3 

q l = P 1 P 1 ' 

q 2 = B 2 P 2 ' an< ^ i 

q 3 = 6 3 P 3 ' 

where P . , q . and 8 . are defined 

i tr i 
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Appendix A 

Development of Equival ent Spring Constant 
The net lift force acting on the balloon is 

Fl = f H k $ ( ff, ■ 0 . 

If (a) the mission takes place in the isothermal 
region of the atmosphere, and 
(b) the helium undergoes an adiabatic process; 
then p n 

.ml?- 

Po 


{A- 1) 


fa = ha. (e ° ) ; 




A -U (e R ) 




(/I -2) 


Substitution of equation (A-2) and (A-3) into equation 
(A-1), gives 


r l " 


Aa 


o 

j OH 


/ 


(A-q 


The equivalent spring constant is obtained by 
differentiating equation ( A— 4 ) with respect to Z, then 
setting Z=0 and considering % to be a constant, 


i.e. 


K - _ PLB 

ci?_ 


L-\ _ 




Fn Vtt $ Pc ( i _ _L \ 
PoH Pc V *« J 


(A-5) 


2=0 
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Appendix B 

Development of Viscous Brag Torques due to 
Rotation of a Spherical Balloon 
(a) Velocity Distribution of Fluid Across the 
Boundary Layer Thickness 

assume (1) parabolic velocity distribution( see figure ) 

(2) a, b,c are determined by the boundary conditions, 
V f = r 0 UJ Sin e at y=0 

- 0 


a 

*2 ” 




at y=£ 
at y=£ 


2 HU) 5 >r>0 

I 


z 


H up $i!2 i 

S' 1 


f. (B-0 


(t) Shear Stress on the Surface of the Balloon 

'ii (ri)r=r 0 , S<<fi 

= -/«• 0"Jr ( ))r=io 3 -4^ r -r,-7“ $$) r „ (S i'D 


( 1 

£ 

14- 

\ CSit 1 6 


0 

( r -2-( 

Jfe 


V <5>f 1 

r 

J/ 


Substitution equation (B-1 ) into (B-2) yield: 




_ -Stfi d 

6 


( 9 - 3 ) 


(c) Viscous Drag Torque 

fTT 

l 

Substitution for ( r r ^) ( _ , we obtain 


At* j (r^) rsri f2Trr;S;neJ(r e 5;/.6;C»7 **) 
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Appendix C 

Development of Inertia Drag Tbrque due to 
Unsteady Rotation of A Spherical Balloon 
in A Viscous Fluid Media 


Average momentum of fluid per unit mass (see figure 9) 

= j \°V f dj si.is 

Angular momentum (11^) of the fluid 


H = \ (average momentum per unit mass) (arm) dm 

~ Y k u).5>('aQ (r 0 St* o ) f (ztt r 0 su\ o') *bcU {As ,.u V^ioY 1 " 

Jo ' / 

12. 5 b (fa V** \tf u To 4 * 


Inertia Drag Torque (N^) is given as 

U T - YY = (f uV^rY 

d-t 
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Appendix D 

■Velocity Expressions for Subs ystems 


Mass 

Vt '= ( x - r, cif, > 4, ) i -i 

( j rr l c(pc((' l ) 6 r qM). )} + 

( j +r,cwhuy e, + r, $($)<: '<*,)% ) | 


Mass m 


7 


v. 


( x -^< 1 ^ 10 % ) i 
4 f y * ifictfiMt ) 4 -gfi sifjmrft 


Mass m. 


V 3 


( j - X r;c(y £ ;y ; , ) i 

s3 


t 


( j + £ »} c ) c ^ - z r i sips ( o L ) / . ) j 
( J i /, r t c (Vj 3 fa; f v 4 L QS 4i ) J- 



( T>'i) 


(' 0 -?. ) 


(t> " 3 ) 



iH I III IIIINIIII 
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Appendix E 

Expressions for Kinetic Energy of the Balloon and 
Subsystems 

(a) For Balloon 

jC _ TWo ( x*+ Lj\ j* ) f ~ ( ) 

t ~T ( f, J-2'SCf',) 1\ j (£-1) 


(b) For 

T< ,j,* + y f ^#*4 ppp,) tf-mcty) x f, 
t inaf'jcje,) ^r^r r 5(y f ).s(e l )^f, + -2nc(i;)S(^$ s, 

+2 qsifxcfe,)^ <j> J 4 ^S 2 (p + & + *$($) Gt f,) (E-2) 



(c) For 


r 




• 


4 




4 i r> 
1 / 


,J y< t 2r}ccff))zs(ps(e r e^) e,% 


f 2 n Sip r 2 5 (%)c(e r 6i) V,t+ 2 nc($)£<xfjc(e r %)e& 



+ Jr > c( ft) G cfp f, % i-2 qc('l l )c(c l )j-h -2rj6(pJ (o,)^f, 
i*bc(%)c(e x )2 i -2 t^(%ysp) j ^ $, 

i *r, s(pc[e ( ) j f, +x Gc(f 2 js(e x)i j % 

t ( PP( p 4 ff+ 2 f-, 4^(%) ) ~ 3 ) 
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(d) For 

T' = x ( X 2 f y "i V * tjccfi > /l-5( y; ) -5" fe, - (±) $1 ( ii 

1 2 1 7 iS- 5 % j c> (tf r ej Vj i f 2 r, C( ft ) 5 Oft , } " h -a ; 4 4 
-;> r r s(ft> i;cty)S(Q r *j>) i, A i\L r ic i( p4 'J*X/} C ( %) 


3 

y S 1' 

3 


c l%)-? 



■+*nrty l )tzc(fr)+ l i + x fa.^. 

+ ■? j J £ 4 ; ~t a j 2 r- 

+ 2^ %^p£ t '< 4 r ^ ^ J ^ ^ 'j ^ )C fe"~®3 ) 

t> r 3 V^/fj ki C( P ^(W'^3 $S C (%) Ojj 

t 2*3% C H) ZQ% C (%) t ^r(* 3 s 2 ( 'fs ) f 4 " 

+ a %•>(£) fj ,) 


(£- 4 ; 
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Appendix F 

Development of Expressions for Generalized Momenta 


?x* Jo” 1 ** * jpWW * >% - . 

Pj- 4**-7 4 

+ >Vi £(%)c(e } )h- '4, 

- »"i r i s(%)5(t> 3 )% 

, a £ 

U - i- Z^qcOf^sfe,) {\j t CcH('J Jfa) Z 

4 >*jr 3 c(y 3 )j(fc\) 4 t Z^iqJl i 

+ v\-ir v 5 (^)c(t> 3 ) % 


( f-t ) 
( F- 2 ) 


(F~ 3) 



(iv.icZf,) i i«sY?,; t^r^ZOf,) ) 0 f 

< - I 

+ £ ? Z r < c %> ( Z^)c(e r oZ ^ c{i h)Z(h) c (z-h) 

5. 4 - 

3 . . ; 

+ + *413 5Y ii)t]c(^)S(e r e f )^ 


(F-4) 


Pt 2 : = ( im. r 2 c(f x )r(e a )j •fZW ( GC(^)o^jj 

i X 3 r 3 c(^ ) r,o(l)C(e } -e,; 4 - A»y 17 J<y,;G ft 

•* —I? 

f >»*3^-^‘fe;r i «-’^;«s{’v03W.5 + 3*3 ^ ^y-O 


(F -5) 


R> 3 - »<$%«%)?(*!) j -f 4 Z H 3 r 3 a<f ; )n^/i) c fo-%) 4 (^"4 

4 ( i 3 , 4 ^) 4 Mj n c ? ( 44 ) ^ 4 1.33 S(% ) f‘, 

';| >)1 4 , 3 ' °3 ) t 



(F-7) 




i - 1 

' x -z^qsipsc c,)ij 4 .f^n^yW'Vv ^ 

t (t.,- 4- iw,n 7 )+, 4 £c(£) f J»V/ Sty )G5(%)c(^)) % 

+ ( >*3 C» -^7 ; 7 7(7, k: <4,- c_ ? ) i »< 3 r ? c ff; )qri ^) ) % 

■1 

■ 4^' nstyJG c(^)S(9 r e^) 6, - >« 3 1] % 


'X 3 h ^v(f)x i^-izsn rf-g^) 

i J; 4V n ay,- ; £*5 ( Yi -e> j 4-, - >" 3 >S < ty > 70 ; (y ; K ( y - e> ) <7 

+ ( J>4 qS(f () qjOq ) c(p j -e x ) 4 jt ^r7*fy£C(fj) ft 
4 .§ H 4 V; 4 (»t } )}jcy, jqs(%)c(q r e 3 ) 4- m 3 ^ ) g e (V o) Vs 


p^3 - - >< 3 >3 > x - m 3 i $s (%$(%)$. rr ?; 

4^ 3 f3 J( ( /3) i7c(ypj(e | -e :3/ ) 6, 4 ^Tr^j^c ^ 

4 0 3 15 ~S(f s ) \f s(f, K ( 4 - fr* ; + >*13 5 <-'44 ’ 4 c( 'h > ) fi 
4 ( >!( 5 (3 5 (% ) 17 5f%) c<t\- e 2 ) t w* 4 C(% u; c(% j)^ 

+ w j 4 +3 • 

pf, = (ic 5 s(‘/ i )4i, 5 5(^.>; 4 -#- ( 1*3 f I B H, } ) 4, fF-/o; 

4 - 4 ?3 j:( 4) e 2 


P<f>3 — 1.33 O3 4 c p 
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Appendix G 

Development of Ex pressions for State V a riable 
Equations 


fit =& V'h 4 - V j\ V H ( x - V, ) - Kd /A/u [ (i }f ( j - £) + y j ¥(x~v< ) 

P^ s pn 4 - i/fl. v h ( - vTy- ) 


/C. £ - i 


T> 


jj z u-v) 

{P-/ 2 ) 

-,-J- 
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y, t £»^n^f,)r;3(%) c (*/-&>) * t l k 
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- h, 3 g o ^yj )3( e r e s > % % - c ^h > i c (e i J 4 ^ 

' hh i Y ^'( i f^ r \ r - t 'i,) r ( fx l- ( '^) H, f 3 J + wx*'(* t ) 


' K, r 


I '0 


( + *f ,*f >i)F [ (fn/^ ) j ^ (/„ /v, >’ J ( 4 +- j) 3(f, > ) 

- K. [y,/v.fif (/h/'hF] rF( h'+Hi±tf)*(*t^rh c(, f<> f $>* ( V) 

(V'MW 




+ f - — tt P r 5 


3 ~t Jk t; 


-i \ + |f/\, mV' 

/ r - 1 (v> G-. V *• (V> C,.^ / 

ttJ h rf<(\, r t ) 


t 

/ 1 * * 

1 1 2 (kr 0 t o(: v : r ( ; J 

( 0/f HAciw-t-h-stf,)) 


4 
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Appendix H 

Expressions for the Matrices , C.. and 
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Expressions for the Matrices 
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Expressions for the Matrices 
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Appendix L 

Expressions for the Matrices Y< , S/ 
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Appendix M 
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